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Most of our day-to-day experiences are mediated by 
light and charged particles, and in particular its interac¬ 
tion with electrons. To the best of our knowledge, the 
physics of a system of photons and electrons is indepen¬ 
dent of the velocity and orientation of that system in ab¬ 
solute space, nor is it locally dependent upon where that 
system lies in a gravitational potential. These symme¬ 
tries, respectively described as local Lorentz invariance 
(LLI) and Einstein’s equivalence principle (EEP), are 
fundamental to our modern understanding of the stan¬ 
dard model and general relativity. It is possible, however, 
that these symmetries are not exact at experimentally ac¬ 
cessible energy scales, thanks to spontaneous symmetry 
breaking or other physics at high energy scales [ 2 , 0 j. 
This possibility has driven many experimental tests of 
LLI and EEP 0 ], and motivated the development of 
phenomenological frameworks that can quantitatively de¬ 
scribe the effect of LLI- and EEP-violation on known par¬ 
ticles and fields. One such framework is provided by the 
standard model extension (SME) 0,1 , which has been 
used to analyze a wide range of experiments [@]. The 
SME augments the standard model Lagrangian with all 
combinations of known particles and fields that are not 
invariant under Lorentz transformations, but which pre¬ 
serve gauge invariance, energy and translational invari¬ 
ance, and the invariance of the total action 01 - These 
terms are parameterized by Lorentz tensors that are col¬ 
lectively known as LLI- and EEP-violating coefficients, 
and are further subdivided into ’sectors’ that deal with 
terms involving a particular particle. In this Letter, we 
focus on tests of the electron-sector c M „ tensor, which 
modifies the inertial energy of electrons according to their 
direction of motion. 

Spectroscopy of neutral dysprosium atoms has already 
led to one of the world’s most sensitive tests of elec¬ 
tronic LLI and EEP 0|. More recently, a still more sen¬ 
sitive measurement of the electron ccoefficients was 
obtained by engineering the quantum state of a pair of 
trapped Ca + ions Q, extending precision tests of elec¬ 
tronic LLI past the electroweak (relative to the Planck 


mass) scale. Both of these experiments operate at or near 
the interrogation-time-of-flight or atom (ion) shot-noise 
limit. In this Letter, we consider the possibility of using 
rare earth ions doped in a crystalline lattice to perform 
similar measurements of the electronic c p V . Rare-earth 
ion-doped crystals offer substantially larger ion-number 
densities than are available in atomic gases, and interro¬ 
gation times comparable and potentially longer than are 
possible in ion traps. The relevant 4f orbitals are well 
screened from one another and from the fluctuations of 
external fields, yielding the sharp, stable optical tran¬ 
sitions that rare earth ion-doped materials are known 
for at low temperatures. The strong crystal field pro¬ 
duces optical-frequency splitting of the otherwise degen¬ 
erate free-ion \J,M) states. This splitting is far larger 
than could be produced by an externally generated mag¬ 
netic or electric field, and as a result, is very stable with 
respect to external field fluctuations. These properties 
make rare-earth ion-doped crystals extremely advanta¬ 
geous for fundamental symmetry tests, complementing 
recent tests focusing on the resonant modes of bulk crys¬ 
tals @. Our work raises the novel prospect of using solid 
state systems to test LLI with greater precision than can 
be achieved by spectroscopy of free particles. 


I. THEORY 

In the SME, spin-independent violations of LLI for 
electrons generate a linearized perturbation of the elec¬ 
tron v tensor to the Dirac Hamiltonian which may be 
written in natural units as mmm 

5H = -{c 00 Sj k + c jk )a 3 p k - ( c 0j + c^p 3 - c 00 ml3, (1) 

where is a symmetric, traceless, constant background 
tensor, a 3 = "/ 0r y 3 and j3 = 7 0 are the usual Dirac ma¬ 
trices and Sik is the Kronecker symbol. In general, one 
can define a coordinate transformation which maps the 
elements of c p V to zero, at the cost of generating new LLI- 
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violating terms in the photon (and other matter) Hamil¬ 
tonians. Such terms can in turn be constrained by mod¬ 
ern Michelson-Morley experiments HI. Since our choice 
of coordinates is arbitrary, it follows that measurements 
of Cfu, are equivalent to Michelson-Morley tests, as both 
are only sensitive to the differences in LLI-violation in the 
photon and matter sectors. At first order, the last two 
terms of Eq. © do not contribute to shifts in the tran¬ 
sition energies between different electronic bound states, 
and may be omitted. We can therefore write Eq. © in 
terms of the spherical components of irreducible tensor 
operators as: 

SH = -C^T^ - Y (2) 

q=-2 

In the spherical operator form used in this Letter, the 
elements of c MJ/ are written as [H, 13 

r>(°) — ~ i 2 „ _2U„ r (2) _ l/o_ _ „ i 

Eg — Coo + 3 Cjj 3^" Coo Lg — g [OC ZZ Cjj) 

C±l = i( Czx i ic z y ) L f g_2 = 2^ xx Cyy A lic x y) 

where doubled roman indexes indicate a sum over the 
spatial components of c^^, and U is the local Newtonian 
gravitational potential. 

The problem of a rare-earth ion in a crystal field has 
been well studied fl5l4l7l |- The total Hamiltonian of a 
rare-earth ion in a crystal field may be written as H = 
Hf + V, where Hf is the free ion Hamiltonian and V is 
the electrostatic potential due to the crystal environment. 
The states of free ions are spherically symmetric, and 
designated by their total angular momentum J and its 
projection M. When the ions are inserted into a crystal, 
the ambient crystal field V breaks spherical symmetry, 
partially lifting thej2J + l)-fold degeneracy of the free 
ions’ energy levels 3- The crystal potential V can 
be written as 

( 3 ) 

k q——k 

where = y / 47r/(2fc+ 1 )Y^ q are normalized spherical 
harmonics, and Yk q are spherical tensor operators of rank 
k. The summation over k is restricted to even numbers 
because the contribution of V to leading order energy 
shifts must come from its even-parity components, and 
k < 6 due to the triangle condition for spherical harmonic 
integrals (since l = 3 for rare-earth ions with configura¬ 
tion 4 f N ). The number of terms in Eq. © may be further 
reduced using the discrete point symmetry of the crys¬ 
tal. The B^ coefficients, also known as the crystal field 
parameters, depend on the structure of the crystal and 
the electronic wave functions’ radial components, and are 
determined by a least-squares fit to the experimental en¬ 
ergy levels of the ion in the crystal. The crystal field 
potential is assumed to act only on the electrons in an 
open shell, i.e. 4 f N for rare-earth ions. 


We obtain the crystal field-induced energy split¬ 
ting from the secular determinant \(JM\V\J'M') — 
X6jj>6mm' | = 0 acting on the free ion states | JM). Diag- 
onalization of V separately within each J-manifold yields 
the split eigenstates |i/>), so that \ip) = um\JM). 
Thus the non-zero matrix elements of m are restricted 
to (JM\Tq\JM'), with q = M — M' , and |g| < 2. 

The LLI-violating correction 8uj nm = (SE n — SE m )/h 
to the transition frequency u> nrn between each pair of lev¬ 
els n and m of the ion in the crystal field is a linear 
combination of the spatial components of the c ^ tensor, 
which is itself a frame-dependent quantity. This energy 
shift also varies as a function of the ion’s position in an 
external gravitational field, although as we will see, this 
effect is smaller than the frame-dependent phenomena for 
the transitions of interest. To uniquely specify the value 
of Cy ,we must also specify the inertial frame in which 
it is defined. This frame is typically taken to be approxi¬ 
mated by the rest frame of the Sun: specifically the sun- 
centered celestial equatorial frame (SCCEF), denoted by 
coordinates (T, X, Y, Z ), while the local laboratory frame 
coordinates are denoted as ( t,x,y,z ). For a terrestrial 
laboratory, the lab-frame values of the tensor’s dominant 
spatial components c*^ b depend upon the orientation of 
the lab with respect to the SCCEF, and thus modulate 
with characteristic frequency ~ 2 x 2-7r/(23 h 56 min), 
or twice every sidereal day. The value of the anomalous 
tensor in the lab-frame can be related to that in the SC¬ 
CEF via cjf b = A^A ^c^ CEF , where Ais the stan¬ 
dard Lorentz boost plus rotation from the SCCEF to the 
lab frame [7], [l(| Gil • Thanks to the Earth’s orbital veloc¬ 
ity, the boost A^“ mixes the time and spatial components 
of C SCC EF i n to the spatial components of cjf b . This gives 
a measurement that searches for yearly modulations of 
SH access to the parity-odd c® j CEF and the isotropic 
cfCC EF components of the anomalous tensor, albeit with 
a sensitivity that is suppressed by one and two factors 
of the Earth’s orbital boost velocity /3<£ ~ 1 x 10” 4 . 
In what follows, we will focus on the laboratory-frame 
values of c ^ and the corresponding spherical operator 
elements and Cq 2 \ and drop the frame-identifying 
superscript. 


II. RESULTS AND DISCUSSION 
A. One valence electron: Ce 3+ ion in LaCL 

Given the available eigenstates of the rare-earth ions in 
the crystal field, we can easily calculate the perturbation 
Suj nm due to Eq. ©. For trivalent rare earth ions ( R 3+ ) 
in the LaClg lattice (see, e.g., 3: Pg- 149), the crystal 
field has the point symmetry C 3 h , and is determined by 
four crystal field parameters B : 

V = B 2 0 C {2) + B^ 4) + B%C^ + B%(cf + Ci 6 >). (4) 
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The simplest rare earth ion to which Eq. applies is 
Ce 3+ , with configuration 4/ 1 . Following the labeling and 
methods of na , the eigenstates and corresponding ener¬ 
gies of Ce 3+ are presented in Table |T| In the first column 
of the table, labels of the states are taken from the Ref. 
G3- The eigenstates are also distinguished by their crys¬ 
tal quantum number [i. 

The states of the single valence electron 4 f 1 of Ce 3+ 
are linear combination of Dirac spinors with mixing of 
free ion states with angular momenta J = 7/2 and 
J = 5/2. These levels’ LLI-violating energy shifts ([2]) 

follow from the Wigner-Eckart theorem and the reduced 

( 2 ) 

matrix element of the tensor C q . In terms of the expan¬ 
sion coefficients and radial integrals I(k' , k) given in the 
Appendix, the states’ shifts are linear combinations of 
radial integrals /(3,3) = —50.0, /(—4, —4) = —49.33, 
and /(—4,3) = —49.58 in atomic units (a.u.). Here 
the radial integrals can be obtained from the formu¬ 
lae I 5 and I\ in the Appendix, and are taken over 
Hartree-Fock (HF) wave functions. The total shift of 
each state is presented in the rightmost column of Ta¬ 
ble HI The largest relative energy shift is that between 
the ground 1 1) state and the low-lying | II) state, with 
Sojiji = (2tt)Cq 2 ' 1 (2.76 x 10 16 Hz). A similarly large LLI- 
violating energy shift is observed for the 1 1 ) to | b) tran¬ 
sition. Details of this calculation are presented in the 
Appendix. 

In contrast to the case of neutral dysprosium, the 
contribution of the scalar component of the LLI- 
violating perturbation in Eq. © to the ions’ transi¬ 
tion energies is smaller than that of the tensor opera¬ 
tor T 9 (2) Q. This occurs because the low-lying ion exci¬ 
tations are largely between states with the same quan- 

I 

6 E{J,ri = -C™<jiJ || T® || M J) 


where clm are the coefficients for the wave functions pre¬ 
sented in Tab. El Inspection of Tab. nn reveals that the 
maximal change in the ions’ transition frequencies due to 
LLI-violation is 5ub,c = 2 ttCq 2> (5.6 x 10 16 Hz) for Nd 3+ 
and <5w Cj( . = 2ttCq 2 \i1 x 10 16 Hz) for Er 3+ . As in the case 

of Ce 3+ , the scalar shift proportional to Cq 0 ' 1 is expected 
to be comparatively negligible. 

A dedicated experiment measuring the THz-scale en¬ 
ergy splitting between the |c) and |e) states of Er 3+ at 
the level of 1 mHz would be sensitive to Cq as small 
as 1CT 20 . To estimate the reach of existing experimental 
measurements of Er 3+ transitions, we have also consid¬ 
ered optical transitions in Er 3+ :Y 2 SiC >5 [2C|. Because 
of the lower symmetry Cf ft of the crystal Y 2 Si0 5 field, 
each J manifold is split into J + 1/2 doubly degener¬ 
ate states (2l|. We focus specifically on the Z\ —\\ 


turn number n and total angular momentum J. In the 
non-relativistic limit, the scalar operator Tq is propor¬ 
tional to the sum of the bound electrons’ kinetic energy 
E, pf/2m. Using the virial theorem in a Coulomb po¬ 
tential, the electrons’ binding energy is approximately 
equal to their kinetic energy. The LLI-violating change 
in the transition energy between bound states with the 
same n and different J is therefore expected to scale as 
Sujnm ~ Uq°^( 6 x 10 13 Hz). This, combined with the com¬ 
paratively small range of variation in U/c 2 accessible to a 
terrestrial laboratory (~ 10 -10 in the Sun’s potential over 
a year), implies this ion could serve as a stable reference 
standard to compare against a more sensitive transition 
(such as those offered by dysprosium 0 ) in a null-redshift 
test of EEP. 


B. Several valence electrons: Nd 3+ and Er 3+ ions 
in LaCR. 

We can perform a similar calculation for rare earth 
ions with several valence electrons in LaCU- Using the 
energy levels and approximate wave functions for Nd 3+ 
and Er 3+ available in the literature {n|, we find that 
these ions have respective ground state configurations 
(4/ 3 , 4 * * * * * * /g/ 2 ) and (4/ 11 , 4 /i 5 / 2 ). In this case, the rele¬ 
vant electronic wave functions are linear combinations of 
Slater determinants of HF orbitals. As before, the fine 
structure manifolds of these ions are split by a crystal 
field with point symmetry C^hi and as before, the con¬ 
tribution of the Cq 0 -* anomaly to the observed transition 
frequencies is much smaller than that of Cq . Repeating 
the analysis of part III Al we may write 


3E M a 2 M M 2 * - J(J + 1) 

i/(2 J + 3)(J + 1)(2 J + 1) J(2J - 1) 


(5) 


transition between the lowest energy levels in the ground 

4 -^15/2 and excited 4 /i 3 / 2 manifolds [20J. Using Eq. ©, 

the LLI-violating perturbation proportional from Cq ' to 

this transition’s energy can be obtained from the reduced 

matrix elements: {/xJ || T^ || / iJ) = —67.097 a.u. for 

J = 15/2 and —56.188 a.u. for J = 13/2. As for the 

case of the lowest level of Er 3+ :LaCl 3 , we have taken 

the weighted sum E m a M-^ 2 ~ 10 f° r both levels Z\ 

and Yi (see Table El- This yields the frequency shift 

1 JwzjVi = Cq 2 " 1 (1.03 x 10 16 Hz). Though not considered 

here, we note the C// component of the LLI-violating 

tensor might also contribute to buxziY i- 
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TABLE I. Experimental crystal-field splittings, calculated wave functions, and LLI-violating energy shifts SE for Ce 3+ (4/ 1 ) 
ions in LaCR. Data from reference [13 ■ Bg = 129, Bg = —329, Bg = —997, and B f = 403 in cm' 1 are used in calculating 
Eq. ©. Note that the eigenstates are doubly degenerate in non-magnetic crystals. 


Config. 

State 


Wave functions \ijjy) 

E (cm 4 ) 

SE {C^ 2) x 10 6 cm 1 ) 

^F 7 /2 

1 d) 

±3/2 

0.99924|7/2, =p3/2) ± 0.03905|5/2, q=3/2) 

2399.5 

1.54 


\c) 

±1/2 

±0.82174|7/2, ±7/2) ± 0.56692|7/2, =p5/2) ± 0.05785|5/2, =p5/2> 

2282.6 

-2.52 


|6> 

±5/2 

=F0.99446|7/2, =f 1/2> ± 0.10511 5/2, =pl/2> 

2208.6 

2.47 


a) 

±1/2 

±0.56659|7/2, ±7/2) =f 0.82356 7/2, =p5/2) ± 0.02709|5/2, =p5/2) 

2166.0 

-1.56 

"F 5/ 2 

1 HI) 

±3/2 

0.03905|7/2, =p3/2) =p 0.99924|5/2, =p3/2) 

110.0 

0.62 


1 II) 

±5/2 

±0.10511|7/2, =pl/2) ± 0.9944615/2, =pl/2> 

37.5 

2.61 


1 1) 

±1/2 

=F0.06298|7/2, ±7/2) =f 0.01072|7/2, =p5/2) ± 0.99796|5/2, =p5/2) 

0.0 

-3.16 


TABLE II. Experimental crystal-field splittings, calculated wave functions, and LLI-violating energy shifts SE for the ground 
levels Nd 3+ (4/ 3 , 4 /g/ 2 ) and Er 3+ (4/ 11 , 4 /i 5 / 2 ) in LaCL. Data from Ref. [Ilf], in terms of crystal quantum numbers y from 
Ref. ED . 


Species 

State 


Wave functions | tp n ) 

E (cm” 1 ) 

SE {C ( 0 2) x 10 6 cm 1 ) 

"NcF 7 

l e > 

±3/2 

0.558| ± 9/2) - 0.830| =p 3/2) 

249.4 

0.05 


Id) 

±5/2 

0.936 ±5/2) -0.351 ±7/2) 

244.4 

0.16 


c > 

±3/2 

0.830 ±9/2) ±0.558 ±3/2) 

123.2 

-0.83 


1 b) 

±1/2 

1 ± 1/2) 

115.4 

1.03 


|a) 

±5/2 

0.3511 ±5/2) ± 0.936| ±7/2) 

0.0 

-0.42 

IF' 7 

\h) 

±1/2 

0.905| ± 13/2) ± 0.356| ± 1/2) ± 0.232| =f 11/2) 

229.31 

1.31 


1 9) 

±1/2 

0.820 ± 11/2) ± 0.429 ± 1/2) - 0.379 13/2) 

181.04 

0.45 


\f) 

±3/2 

0.116| ± 15/2) ± 0.764 ± 3/2) ± 0.635| 9/2) 

141.61 

-0.96 


|e) 

±5/2 

0.662| ± 7/2) ± 0.750| ±5/2) 

113.7 

-1.08 


Id) 

±1/2 

0.192| ± 13/2) - 0.830| ± 1/2) ± 0.523| ± 11/2) 

96.52 

-0.26 


c) 

±3/2 

0.925 ± 15/2) ± 0.150| ± 3/2) - 0.349| =f 9/2) 

64.27 

2.58 


|6> 

±3/2 

0.362| ± 15/2) - 0.628| ± 3/2) ± 0.689| =F 9/2) 

37.91 

-0.29 


|a) 

±5/2 

0.750 ± 7/2) — 0.662| ±5/2) 

0.0 

-1.02 


C. Magnetic field effect 

Stray magnetic fields can produce major systematic 
errors in Lorentz symmetry tests. We have therefore es¬ 
timated the effects of magnetic fields on such tests which 
use Ce 3+ and Nd 3+ in LaCl 3 crystals. Detailed stud¬ 
ies of the Zeeman effect in Ce 3+ ions may be found in 
reference [22J, wherein the first order Zeeman splitting 
of doubly degenerate states may be found. Though the 
ground state degeneracy is lifted by the 2.1 /tb (where 
/iB = 0.467cm _1 /T is the Bohr magneton), the mag¬ 
netic interaction between states in different doublets is 
restricted due to the crystal symmetry and the selection 
rules for Ml transitions. For Ce 3+ , the quadratic shift 
has been estimated using the dominant components of 
the ground state |/) and the low lying state | III) at 110 
cm” 1 (see Tab. [I]). Since the magnetic quantum numbers 
differ for these states, transitions are only possible by way 
of the x and y components of the magnetic moment op¬ 
erator. Contributions from other levels corresponding to 
the 2 F 7 / 2 term are suppressed due to the larger energy 
splittings. Thus we estimate the second order shift to 
be 0.918/(110cm _1 )/x|. This corresponds to a quadratic 
shift of 5.46 x 10 7 Hz/T 2 for the ground state. The same 
calculation can be done for the | II) state, for which we 


obtain 13.2 x 10 7 Hz/T 2 . 

The Nd 3+ ion has three valence electrons in the / shell, 
and so we limit ourselves to a rough estimate, which may 
also be applied to other rare-earth ions. The total an¬ 
gular momentum J = 9/2 is bigger than the total spin 
S = 3/2, so the magnetic moment of the ion is domi¬ 
nated by the angular momentum operator J. Hence the 
matrix elements for transitions within a given multiplet 
are easily obtained. As a result, the quadratic shifts are 
expected to be ~ 30 x 10 7 Hz/T 2 for the ground level |a) 
and ~ 18 x 10 8 Hz/T 2 for the first excited doublet |6) 
(see Tab. ITT1). 

Using these estimates, we may determine the extent 
to which magnetic fields must be controlled to suppress 
their effects to below that of an LLI-violating Cq with 
order 10 _2 °. For Ce 3+ ions, the magnetic field must be 
stabilized to be below 4.7 x 10” 6 T, or 47 mG, while 
for Nd 3+ ions, a field of no more than 15 mG should be 
sufficient. We further note that larger magnetic field fluc¬ 
tuations, and larger DC fields, are in principle tolerable 
for a test of LLI, so long as they fluctuate on timescales 
that are sufficiently different from the modulation peri¬ 
ods of the laboratory’s orientation and boost relative to 
a fixed inertial frame ( e.g. as the SCCEF approximates). 
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III. CONCLUSION 

We have demonstrated that solid state systems, and 
particularly the ground state spectrum of rare earth 
ions doped in a crystalline host can be used to per¬ 
form sensitive tests of LLI. We have taken advantage 
of existing work on the spectrum of Ce 3+ pa, Nd 3+ , 
and Er 3+ [23|, HH to perform an explicit calculation of 
these ions’ energy shifts in response to LLI-violation 
when doped in LaCl 3 . The energy levels and sensitiv¬ 
ities of Ce 3+ , Nd 3+ , and Er 3+ ions are expected to be 
similar when doped in different crystalline media. Er 3+ 
is a particularly interesting case, as the optical coher¬ 
ence of the J = 15/2 —y J = 13/2 transition is par¬ 
ticularly long-lived at 4.4 ms for 0.001% doping concen¬ 
tration Er 3 + :Y 2 Si 05 [25|. An experiment that is sen¬ 
sitive to a 1 mHz orientation-dependent modulation of 
the Z\ —> Yi optical transition in Er 3 + :Y 2 SiC >5 could 
measure spatial components of the electronic c ^ ten¬ 
sor as small as 10 ” 19 , improving upon existing limits by 
an order of magnitude [7|, [ 8 |] . An experiment measuring 
orientation-dependent modulations of the THz-scale en¬ 
ergy difference between the |c) and \e) states (see Tab. ED 
at the mHz would be more sensitive still, probing the 
LLI-violating c ^ tensor at the level of 10~ 2 °. We note 
that dynamic decoupling techniques [26| . which switch 
between states with quantum numbers of equal magni¬ 
tude and opposite sign do not suppress the LLI-violating 

( 2 ) 

signal proportional to the quadrupole component of Cq ' 
of C/j, u . Other rare earth ion-doped materials may also 
prove to be useful for testing this and other aspects of 
LLI, and are a promising area for future work. 
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APPENDIX 

In what follows, we restrict the crystal quantum num¬ 
ber to — fi because both degenerate states get the same 
shift from the LLI-violating perturbation. The eigen¬ 
states |7), |a), and |c) are presented in Ref. 0 : 

I^-m> = £17/2, -7/2} + r?|7/2,5/2) + C|5/2,5/2), (6) 

with expansion coefficients £, rj and £. Similarly, the 
remaining eigenstates | II), \III), | 6 ), and | d) are: 

IlM =£|7/2,m) + C|5/2,m), (7) 


where m = 1/2 or 3/2. These basis states are approxi¬ 
mated by the relativistic, four-component spinor Hartree- 
Fock (HF) orbitals of the free Ce 3+ ion. 


_ 1 / fn K (r)Q. K m{6,(t>) \ 

r \iag nK (r)Cl — Km {0,4) ) 


( 8 ) 


where the non-relativistic two-component spinor is de¬ 
fined by 


4) 




(9) 


and k = =f(/ + 1 / 2 ) (for j = l ± 1 / 2 ) is the unified quan¬ 
tum number denoting angular momentum and parity. 

For transitions between states with the same quantum 
number n and total angular momentum J, the main con¬ 
tribution of the LLI-violating perturbation comes from 
the tensor Tq ’ operator. Its matrix elements may be 
written as 


{n'K'm'\Tjj 2 ^\nnm) = (—l )' 7 


f 2 j 

—m! q m 
x ( k ' || Cf || k)I(k', k). 


( 10 ) 


where I(k\k) represents a radial integral [27]. The re- 

( 2 ) 

duced matrix element of the tensor Cq ' is 


W II C f II «) = 


(-1 ) j,+1 *VWJ\ 


J J K 

_i I o 
2 2 u 


(ii) 


where \j',j\ = (2 j' + 1)(2 j + 1). The energy shifts due 
to LLI violation can be calculated using equations m, 
and CEB- In terms of the expansion coefficients and radial 
integrals I(k',k), the shifts of the states 1 1), |a), and |c) 
in Eq. © are 

6E = ^-((7e + y 2 )l(-4,-4) + 

+6C 2 /(3,3) + 2v / 6C^(-4,3)). (12) 


Two other states | II) and | 6 ) have the shifts 
( 7 ( 2 ) 04 

SE = -^(5£ 2 /(-4, -4) + —C 2 7(3, 3) - 

-^«/(”4,3)) (13) 

with m = 1/2 in Eq. ([7|l. while the states | III) and | d) 
have the energy shifts of the levels 

C ( 2 ) c 

SE = -^(3£ 2 /(-4, -4) + -C 2 /(3,3) - 

4,3)) (14) 

Here the radial integrals can be obtained from the for¬ 
mulae -Z 5 and I\ below, and are taken over HF wave func¬ 
tions. 
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The radial integrals included in Eq. (U2I13I14[) are from 
Ref. [j], and are summarized in this Appendix for con¬ 
venience. According to the quantum numbers k' and k, 
the radial integrals take one of the following forms: 


I\ = cah 


i>( 


/Q 1 \ dfnK / 0 . Q\ r dg 

\2k 1 jQn'K' o [2k ~\~ o)J n ' K ' 


(2 K — 1 )(k + 1 ) 


~9n 


r\ I V"' 1 ' I / J ri K r\ 

or or 

f (2av + 3 )k f 
' JnK Jn'K' g n 


for k! = —k — 1, 


for k' = k + 2, and 


I 5 = —2 cah 




dfnK 

dr 


t dg nK 

Jti ' k ,' 

or 


+ 




for k' = k. In the nonrelativistic limit the matrix ele¬ 
ment (</>|cadp fc |(/>} becomes p k / m\<f >). This changes 

the radial integrals in Eq. (9). Only the radial integrals 
differ between the relativistic and the nonrelativistic lim¬ 
its, while the angular parts of the matrix elements are the 
same. In the non-relativistic limit, I\ and I 5 reduce to 


I 2 =cahJ dr^(2« - + (2k + 1)/, 


dr 


dr 


(2k-3)k^ f i (2 k + 1)(k - 1) f 

9n' k' JnK 1 Jn'K' 9r 


for k' = —k + 1, 


h = 4 cah j dr( /„ 


for k' = k — 2, 


h = — 4ca?i / dr g n , K . 


dfnK K + 1 

dr r 


gn 1 k' fn 


r A fdfn'K' 

l ir {— 


h = Ti 

I 2 and I 3 become 

f‘QO 

I 2 = n 2 dr 
Jo 


^fn' k' dfnK ! K^K + 1) ^ f 
“77 1 9 Jn'K' Jn 

or r z 


dfn'K' dfnK 2k -l df % 
dr dr ■ h 


n' K' ry 

r " or 


k(k — 2) f f 

2 Jn'K'Jn 


and I 4 takes the form of 


r *2 / 1 / d/n'K' 5/n/i . 2 k + 3 df , 

A = ft / dr I ——---1-/, 

ar ar r 

(/s7 + 3 )(ac + 1) 


n' K' ry 

r ar 


fn'K' fn 
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